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Abstract 



When a single atom is sandwiched in between two electrodes, an atomic 
tunneling device may be realized depending on the distance between the atom 
and the electrodes. We have performed first-principles pesudopotential cal- 
culation in conjunction with a three-dimensional evaluation of the quantum 
scattering matrix for such a device. We predict the quantum conductance and 

resonance transmission properties of the Si and Al atomic tunneling devices. 
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Exploiting the band gap differences of various compound semiconductors such as GaAs 
and AlGaAs, resonance tunneling devices have been fabricated and extensively investigated 
since the original work of Esaki [|lj. A tunneling device may give rise to negative differential 
resistance due to a quantum resonance, leading to useful electronic applications. In this 
work, we propose and theoretically analyze a different kind tunneling device, namely atomic 
scale tunneling systems. Rather than resonant transmitting through the quasi-bound states 
inside a double barrier tunneling structure, in the atomic system a resonance is through 
atomic orbitals. Since atomic orbitals have energy spacings in the range of eV, the quantum 
resonance will not be smeared out by a room temperature. Hence in principle room temper- 
ature quantum resonance tunneling devices can be achieved. Our work is encouraged by the 
measurements using scanning tunneling microscope (STM) on quantum conduction through 
atomic scale tip-substrate systems |2j, where quantized conductance is indeed observed at 
room temperature. 

We propose an atomic tunneling structure as shown schematically in the inset of Fig. 
(3). The essential part of the structure consists of two metallic electrodes sandwiching a 
single atom in between. The whole system could in principle be fabricated on top of an 
insulating substrate by micro-fabricating two electrodes and placing an atom in between 
using the atomic manipulation ability of STM. Fixing the distance between the atom and 
the electrode, d, vacuum barriers can be established on the two sides of the atom, thus a 
double-barrier atomic tunneling (DBAT) device is established |3j . The idea of using vacuum 
barriers to establish a tunnel junction is widely used in STM |3J, but an DBAT structure as 
an atomic quantum functional device is very interesting indeed. The main task of our work, 
to be presented below, is to predict quantum conductance of the DBAT device for various 
values of d, to investigate the formation of a atomic quantum point contact as d is reduced 
thus diminishing the vacuum barriers, and to examine issues which arise when transport is 
mediated by atomic levels. Although we have concentrated on DBAT devices operated on 
a single atom, the idea can be extended to clusters, molecules, or other groups of atoms. 

To capture the atomic degrees of freedom, we have combined the ab initio pseudopotential 

2 



total energy method with a three-dimensional (3D) quantum scattering evaluation of the 
transmission probabilities ||. First, we solve the ground state properties of the atom and the 
two leads by minimizing the Kohn-Sham total energy functional using a plane-wave basis 
set 0. The equilibrium analysis produces the self-consistent effective potential V e ff(r) = 
SU/5p(r) which is seen by all the electrons. Here U[p] is the total self-consistent potential 
energy while p the electron density. Second, we evaluate the scattering matrix of an electron 
traversing the system defined by V e ff, by solving a 3D scattering problem using a transfer 
matrix technique . From the scattering matrix we obtain conductance from the Landauer 
formula ||. 

We have focused on DBAT devices made of Si and Al atoms and used pseudopotentials of 
Refs. 0|lOj for the core, and the parameterization of Ref. [llj] for the exchange-correlation 



term. The leads are modeled by the jellium model. For Si, a jellium lead has a cross 
section area of 7.25 x 7.25(a.w.) 2 , length L = 30.78a.u., and its charge is specified [12|] by 
r s pa 2.0a.u., mimicking metallic leads. The supercell volume used in our plane-wave based 
ah initio calculations is 21.77 x 21.77 x 2(L + d) (a.u.) 3 . For the Al DBAT device, the lead 
has a cross section area of 8.79 x 8.79(a.w.) 2 , length L = 23. 57a. u., charge density specified 
by r s pa 2.07a.u.. The supercell size for the Al DBAT device is 16.67 x 16.67 x 2(L + d) 
(a.u.) 3 . The DBAT setup with square shaped leads has a symmetry of space group D4h. 
We have used an upper energy cutoff of 8 Rydbergs [13[ . 



For a single adatom on top of a high density jellium substrate, previous calculations 



showed [|14]] that the equilibrium jellium-atom bond length is 2.3a.u. for Si, and 2.6a.u. for 



Al ||12|| . Hence for values of d greater than these values, a vacuum barrier may be established. 



Fig. (1) shows the V e ff for a Si DBAT device with d = 6.9 a.u.. V e ff in the 3D leads is 
essentially a potential well with a depth ~ — 0.50a.u. below the Fermi level of the system. 
From V e ff it is clear that the atom is quite isolated from the leads by the vacuum barriers. 
Reducing d to 3.45 and 4.6a.u., a lower vacuum barrier is obtained as expected. Similar 
behavior is found for Al DBAT devices. In Fig. (1), the sharp peak at the atomic position 
reflects the repulsive atomic core, and surrounding the core there is the usual attractive part 



of the atomic potential. 

Figs. (2,3) show the conductance G for a number of Si and Al DBAT devices which differ 
by the jellium-atom junction distance d. When d takes the value of equilibrium jellium-atom 
bond length, G(E) shows the expected "quantized" conductance |0| for quasi-lD quantum 
wires. However a striking result is the apparent resonance transmission when the atom is 
isolated inside the DBAT device by the barriers. For both Si and Al systems, the larger 
the d, the sharper the resonances. This is because a larger d corresponds to a more isolated 
atom. The resonance peaks must come from the atomic orbitals since the DBAT device 
is operating on the electronic structure of the single atom. For a Si atom, the valence 
configuration is 3s 2 3p 2 . Thus the equilibrium "Fermi" energy for an isolated atom should be 
located at the 3p atomic level. For all the Si DBAT devices, our calculated Ef pa — 0.10a. u., 
which is very close to the higher energy peak position of Fig. (2). This allows us to identify 
the resonance peak near —0.09 ~ — O.lOa.u. as due to the 3p atomic orbital. For the largest 
d we studied, e.g. d = 6.9a.u., the sharp 3p resonance is split into two nearby peaks as 
shown in Fig. (2). It is not difficult to understand this splitting when we realize that the 
3p state is triplely degenerate, and its rotational symmetry is broken by the presence of 
the squared shaped leads which is along the z-direction. Hence for the single 3p z state, the 
resonance is marked by G(E 3pz ) pa 1 in units of 2e 2 /h. On the other hand, the D4h space 
group respects the rotational symmetry in the x — y direction, thus the degenerate 3p x ^ y 
states give G(E 3pxy ) pa 2. The reason that these peaks are slightly higher than 1 and 2 is 
due to the partial overlap of the states, as is clearly seen in Fig. (2). When the lead-atom 
distance d is smaller, such as d = 3.45a.u. and 4.6a.u., the atom is not as isolated. For these 
cases, the resonance width is larger which smears out the 3p splits. As a result, the three 
states of 3p atomic orbital gives a G(E 3p ) = 3 x (2e 2 /h) resonance. 

After understanding the 3p resonance peak, it is clear that the lower energy peak of 
Fig. (2) must result from the 3s atomic state. Indeed, since 3s is a singlet, we have 
G(E 3s ) = 1 x {2e 2 /h). For d = 6.9a.u., the 3s peak is extremely sharp, we reproduce it 
in the inset of Fig. (2). A crucial numerical test of these results is the distance between 
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the peak positions. For a single Si atom, LDA calculation [lfj gives energy level spacing 
between the 3s and 3p states to be 0.245a.u.. Our conductance calculation, as shown in Fig. 
(2), gives the 3s — 3p resonance distance in almost perfect agreement with the LDA spectra, 
with only small differences due to the level splitting and the presence of the leads. Hence it 
is unambiguous that the transport in a single atom DBAT device is mediated by the atomic 
levels. For Al DBAT devices, as shown in Fig. (3), all the results give a physical picture 
which is consistent with that of the Si system, including the quantitative value of the 3s-3p 
resonance peak distance. For different values of d, there is a small but noticeable shift of the 
relative peak positions (see Figs. (2,3)). This is due to the coupling of the atom to the leads, 



and can be understood from previous model calculations |T7|]. Finally, Figs. (2,3) shows the 



formation of a transmissive quantum point contact: by reducing d to the equilibrium bond 
length between the atom and the leads, the resonance transmission crosses over to the usual 
quantized conductance with a quasi-lD nature. 

Certain features of the ab initio results can be understood analytically. First, it can be 
generally shown that for 3D resonance tunneling through a symmetric double barrier 
potential, the peak value of the conductance G(E) for an isolated singlet state is universally 
2e 2 /h. The 3s peak shows this result. Second, using the discussion of Ref. and Fermi's 
golden rule, it's not difficult to obtain the resonance width as 



T = 2vr^ \Mj\ 2 6(E r - Ej) = hJ2 
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where Wj is the tunneling rate from the atomic resonance level E r to the lead level Ej, 
and Mj is the corresponding tunneling matrix element from E r to Ej. Since the tunneling 
rate decreases exponentially with the barrier size, we expect the resonance widths to do the 
same, which is clearly seen in our numerical results. Third, we can understand why the 3p z 
resonance peak width is wider than that of the 3p xy peak. Using a first-order approximation 
for the tunneling matrix pi 



Mj = / x)V lead <p r d\ 
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where Xj is the lead state, r is the atomic state, and Vi ea d is the potential due to the lead. 
By approximating the lead wavef unctions with finite square well wavef unctions, and VJ ea(J as 
a potential well, we find T p jT Pxy = 6.2 for d = 6.9a.u. This agrees well with the value of 6.0 
from our ab initio calculation. We thus conclude, that the 3p z resonance peak width is wider 
than that of the 3p X;V peak, because the overlap of the 3p z state with the lead wavefunctions 
is much larger than the corresponding overlap for the 3p xy states. 

In summary, we have proposed and theoretically investigated atomic scale double barrier 
tunneling devices. The quantum resonance transmission in this system is through atomic 
orbitals, thus conductances obtained characterizes the conduction through a single atom 
and reveals the valence atomic spectra. Finally we comment that the curves of G — G(E) 
are useful for giving a first estimate of the electrical current as a function of a voltage across 
the DBAT device: one obtains the current in the usual fashion by integrating G(E) over 
energy convoluted with a Fermi function. The importance of the atomic DBAT device lies 
in the fact that the quantum resonances are observable at room temperature due to the 
large energy scales (eV) associated with the atomic system. Indeed, with the nonlinear 
conductance behavior of the DBAT device, many applications can be envisioned. 
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FIGURES 

FIG. 1. Effective potential V e ff for a single atom Si DBAT device. The atom-lead distance is 
fixed at d = 6.9a.u.. The vacuum barriers are clearly seen. 

FIG. 2. Conductance G(E) as a function of incoming electron energy E for Si DBAT de- 
vices. Dashed line: for d = 2.3a.u., which is the equilibrium atom-leads distance; dotted line: for 
d = 3.45a.u.; long dashed line: for d = 4.6a.u.; solid line: for d = 6.9a.u.. Inset: resolving the 3s 
resonance peak for the d = 6.9a.u. case. The thin vertical line indicates the calculated Fermi level 
of the system. 

FIG. 3. Conductance G for Al DBAT devices. Dashed line: for d = 2.6a.u.; dotted line: 
d = 3.9a.u.; long dashed line: d = 5.2a.u.; solid line: d = 6.5a.u.. Inset: schematic plot of the 
atomic DBAT device: an atom is sandwiched in between two metallic wires. The thin vertical line 
indicates the calculated Fermi level of the system. 
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